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ABSTRACT 

Weak lensing by large-scale structure allows a direct measure of the dark matter distribution. We 
have used parallel images taken with the Space Telescope Imaging Spectrograph (STIS) on the Hubble 
Space Telescope to measure weak lensing, or cosmic shear. We measure the shapes of 26036 galaxies 
in 1292 STIS fields and measure the shear variance at a scale of 0.51 arcminutes. The charge transfer 
efficiency (CTE) of STIS has degraded over time and introduces a spurious ellipticity into galaxy shapes 
during the readout process. We correct for this effect as a function of signal to noise and CCD position. 
We further show that the detected cosmic shear signal is nearly constant in time over the approximately 
four years of observation. We detect cosmic shear at the 5.1c level, and our measurement of the shear 
variance is consistent with theoretical predictions in a ACDM universe. This provides a measure of the 
normalization of the mass power spectrum <7g = (1.02 ± 0.16)(0.3/fi m ) 46 (0.21/r) 018 . The one-cr error 
includes noise, cosmic variance, systematics and the redshift uncertainty of the source galaxies. This is 
consistent with previous cosmic shear measurements, but tends to favor those with a high value of erg. It 
is also consistent with the recent determination of erg from the Wilkinson Microwave Anisotropy Probe 
(WMAP) experiment. 

Subject headings: cosmology: observations - gravitational lensing - methods- dark matter 



1. INTRODUCTION 

Weak gravitational lensing by large-scale structure has 
become an important tool in understanding the amount 
and distribution of dark matter (see Mellier et al. 2002 & 
Refregier 2003 for reviews). Most cosmic shear measure- 
ments have been performed using wide-field ground-based 
telescopes. However, the measurement of galaxy shapes 
from the ground is limited by atmospheric seeing, while the 
Hubble Space Telescope (HST) has a much smaller Point 
Spread Function (PSF) (< 0.1" as opposed to - 0.8" for 
ground-based telescopes). Several groups have thus used 
HST images to study cosmic shear (Rhodes, Refregier & 
Groth 2000 [RRGI] and 2001 [RRGII] ; Refregier, Rhodes 
& Groth 2002 [RRGIII]; Hammerle et al. 2002). Recent 
measurements of cosmic shear from the ground (Bacon 
et al., 2002, Hoekstra et al. 2002, & van Waerbeke et 
al. 2002; Jarvis et al. 2003; Brown et al. 2003) and from 
space (RRGIII) have been used to derive an estimate of 
the normalization erg of the matter power spectrum with an 
uncertainty comparable to that of more traditional meth- 
ods such as the measurement of cluster abundances (see 
RRGIII for a discussion). 

In this paper, we present the results of a weak lensing 
analysis of images taken with the Space Telescope Imaging 
Spectrograph (STIS) on HST. The STIS data we use have 
been described by Pirzkal et al. (2001). Hammerle, et 
al. (2002) have presented an analysis of a subset (121 
of 1292 fields) of these data, reporting a 4% RMS shear 



at the 1.5cr level. Because there are a number of different 
methods for measuring shear and correcting for systematic 
effects, it is beneficial to compare the results we present 
here to the results of other groups examining data using 
the same instrument. 

This paper is organized as follows. In §2 we describe the 
STIS data used in this study. We detail the STIS parallel 
archive maintained by the STIS Investigation Definition 
Team (IDT) at Goddard Space Flight Center, including 
the methods used to reduce and catalog STIS images. §3 
describes the procedure used to measure the weak lensing 
shear in the STIS data, including the method for removal 
of systematic effects. In §4.1, we present the results, in- 
cluding constraints on cosmological parameters. We draw 
our conclusions and briefly comment on the future of weak 
lensing measurements in §5. 

2. DATA 

The STIS images used in this study were taken primarily 
in parallel mode, meaning that another instrument on HST 
was the prime observing instrument. Thus, the images are 
essentially randomly positioned on the sky, separated from 
the primary observation by 5 to 8 arcminutes depending 
on the primary instrument. This is ideal for a study of 
cosmic shear because many random pointings on the sky 
minimize the error due to cosmic variance. 

There are four stages to the data reduction: (i) input 
image co-alignment for a particular field, (ii) preliminary 
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image reduction that involves bias and dark subtraction, 
hot pixel correction, flat-field division and correction for 
geometric distortion, (iii) cosmic-ray (CR) removal and 
image combination, and (iv) object detection. 

Step (i) involves determining the relative offsets using 
sources identified in dithered images. Copies of the input 
images are pre-cleaned to minimize the effect of CRs on 
the source identification step. Small clusters of pixels with 
high data values relative to surrounding pixels are first 
identified then corrected using interpolation. This pre- 
cleaning may modify or remove real sources, but these 
image copies are used only for the purpose of co-alignment. 
Source Extractor (SExtractor) (Bcrtin & Arnouts 1996) 
is used to identify sources in the pre-cleaned image copies. 
Those sources in common on all images are used to find 
individual image offsets. The average offset is subtracted 
from the individual offsets so that all images will be shifted 
to a common field center. 

For step (ii), we used the Interactive Data Language 
(IDL) version of CALSTIS (Lindler 2003), which was de- 
veloped by the STIS Instrument Definition Team (IDT). 
A one-dimensional fit is made to the CCD overscan region 
and subtracted from each column of the image array to ac- 
count for temporal variations between the observation and 
the two-dimensional bias image. A two-dimensional bias 
frame is then subtracted from the science image. Next, 
the dark image subtraction is performed and hot pixels 
are corrected. The hot pixel tables from the Space Tele- 
scope Science Institute archive list new hot-pixels created 
by the constant on-orbit CR flux since the last dark im- 
age was obtained. For a given observation, the hot pixels 
corrected are those in common between the two lists made 
closest in time before and after the observation. Hot pix- 
els with a dark rate greater than 0.02 counts/second were 
corrected by linear interpolation among the eight adjacent 
pixels. The resultant image is divided by two flat-field im- 
ages: one that accounts for pixel-to-pixel variations, and 
another that corrects the low-frequency spatial variations 
across the field of view due to vignetting. A final hot- 
pixel correction is performed to pre-clean the data before 
cosmic-ray rejection. Using a separate IDL program, small 
clusters of pixels (< 3-pixels wide) with high data values 
relative to surrounding pixels are first identified then cor- 
rected using interpolation. This procedure is especially 
useful in removing hot or bad pixels that would not other- 
wise be eliminated in the CR-rejection step when there is 
no dither. This is the same algorithm described in step (i) 
for removing cosmic-rays before running SExtractor, but 
in this case, the threshold for correction is set to a higher 
value so that valid source pixels arc not modified. The 
correction for detector/optics induced shear (or geometric 
distortion) was performed with bi-lincar interpolation us- 
ing the values given by Malumuth and Bowers (1997). The 
distortion coefficients were derived using the astrometric 
shifts of stellar images. All images are then shifted to a 
common field center using bi-linear interpolation for the 
sub-pixel offsets derived in step (i). 

Once the individual images have been co-aligned they 
are combined with the routine IDL cr_reject .pro (step 
iii). This program emulates the STSDAS task 
stsdas . hst_calib . wf pc . crrej and is equipped to han- 
dle input images with different exposure times. The main 
difference between the two algorithms is in the initial 



step of sky background determination. The procedure 
stsdas . hst_calib .wfpc . crrej uses the modal value of 
all of the image pixels, while crjreject .pro uses an esti- 
mate of modal value of a subset of image pixels. Using a 
subset of image pixels reduces the bias in the pixel distribu- 
tion caused by foreground source flux. The estimate uses 
an algorithm from DA0PH0T I (Stetson, 1987). The aver- 
age and standard deviation of the pixel- value distribution 
are computed, and outliers are removed using an itera- 
tive sigma-clipping method. If the distribution is Gaus- 
sian (uncontaminated by foreground sources) the mean, 
median and mode should be the same. If, after 20 iter- 
ations, the mean and median are the same, this value is 
taken as the sky for the image. If the mean of the distri- 
bution is larger than the median (a non-Gaussian distri- 
bution) the true sky is estimated as 3 x median - 2 x mean. 
The scalar background so derived is subtracted from each 
input image. Pixels are rejected as cosmic-rays if their 
value is greater than an input number of a from a refer- 
ence value. The constant read noise, the statistical noise 
(square-root of the counts), and a noise proportional to 
the counts comprise the a. This last noise component ac- 
counts for differing point-spread-functions (PSFs) or sub- 
pixel image shifts. It is expressed by the "mult_noise" 
variable in crjreject .pro that we set to 0.03. If there 
are 7 or more input images, the first guess at a CR-free 
or reference image is the pixel-by-pixel median of all the 
input images, otherwise the minimum is used. The first 
pixel rejection pass cuts all pixels greater than 8a from 
this reference image. A second pass uses the mean of the 
cleaned images as the reference image and a rejection cri- 
terion of 6cr. The third and final pass uses a 4<r rejection 
criterion on the mean of the cleaned images from the sec- 
ond pass. The sky values that were subtracted from each 
input image before CR-rejection are added back into the 
final set of cleaned images. The final image combination 
is the weighted average of the cleaned images. The weight 
parameters are the sky noise and the read noise. The final 
image pixels are in units of counts per second. 

For step (iv) we use SExtractor (Bcrtin & Arnouts 
1996) for object detection, photometry, star/galaxy sepa- 
ration and to compute positions for each detected object. 
SExtractor subtracts a smoothed two-dimensional back- 
ground with a mesh size of 1" (BACK_SIZE= 20) and filter 
size setting of BACK_FILTERSIZE= 3. The image is then 
convolved with a l'/25 x l'/25 O'.'5-FWHM Gaussian filter. 
The detection and analysis thresholds are both set to 0.7er 
above the measured sky noise. We require that 5 contigu- 
ous pixels meet this threshold criterion to be considered 
a valid source. In order to separate blended sources, we 
set the number of deblending sub-thresholds to 32, and the 
minimum contrast parameter for deblending to 0.005. The 
input stellar seeing disk FWHM, used for star-galaxy sepa- 
ration, is 0'.'07. The photometry reported is taken from the 
MAG.AUTO column of the output catalog. MAG_AUTO 
is the source magnitude measured inside a unique elliptical 
aperture for every object (Kron 1980). The two-element 
input parameter PH0TJUJT0PARAMS controls the elliptical 
apertures. The first element is the k- factor described in 
Bcrtin & Arnouts (1996), and the second is the minimum 
possible radius (in pixels) for an elliptical, or Kron, aper- 
ture. We use the default values of 2.5 and 3.5, respec- 
tively, for the two elements. A 12-pixel annulus outside 



the elliptical aperture is used to determine the local sky 
background around each object. 

We examined 2335 fields imaged in 21 HST programs 
(7781, 7782, 7783, 7908,7910, 7911, 8062, 8064, 8084, 8091, 
8393, 8406, 8470, 8545, 8549, 8562, 8796, 8808, 8870, 
8884, & 9248). Several of these programs are described 
elsewhere: the STIS Parallel Survey (SPS; Gardner et 
al 1998; Teplitz et al. 2003a & 2003b) obtained filterless 
imaging and slitless spectra of random fields on the sky 
and program 8562 (PI Schneider; continued as 9248) was 
conducted to gather data specifically to search for cosmic 
shear (Pirzkal et al. 2001). Of these fields, 1494 pass a 
visual inspection that indicates they are not contaminated 
by stars or large objects. In our measurement of cosmic 
shear, we use 1292 of these fields that pass several other 
cuts described below. 

The selected fields were imaged between August 9, 1997 
and May 15, 2001. We do not use more recent data be- 
cause of possible changes in systematic effects that may 
have occurred when STIS was forced to go to backup elec- 
tronics as a result of a malfunction in July 2001 (Brown 
2001). Each field has between 2 and 86 individual expo- 
sures with total exposure times ranging from 265 to 34,200 
seconds. The distribution of exposure times and number 
of exposures are shown in Figure 1. 
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Fig. 1. — Top: Histogram of the number of individual exposures 
going into each image. Bottom: Histogram of the exposure time for 
each individual image. 

Imaging is done using the 50CCD clear mode, a very 
wide bandpass sensitive to light in the wavelength range 
220 < A < llOOnm and limited by the sensitivity of the 
CCD (Baum ct al. 1996). We co-add exposures that are 
within 5 arcseconds of each other. Only exposures taken 
within 6 months of each other are co-added to produce 
the final image. This limit is chosen to avoid co-adding 
exposures with large differences in their calibration files. 

Using the results of RRG III and the fact that (M) « 
(I) + 1 for typical galaxy colors where M and / are the 
STIS and WFPC2 magnitude respectively, we find that 
the median magnitude M m of the galaxies in each STIS 



field is related to their median redshift z m by 

z m w 0.72 + 0.15 x (M m - 23.0) (1) 

We show a histogram of the median rcdshifts of the se- 
lected fields in Figure 2. 
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Fig. 2. — Distribution of the median redshifts of the 1292 fields we 
use for cosmic shear measurements. The redshift z m is estimated 
from the median magnitude of the selected galaxies in the field via 
Equation (1). 

We visually examined 2335 fields in the STIS archive. 
Our selection criteria were: two or more exposures, un- 
binned pixels (to maintain the full STIS resolution), and 
imaging performed with the clear 50CCD filter. We dis- 
carded fields that have a very bright saturated star, fields 
that are crowded by stars (common for fields with galac- 
tic latitude \b\ < 20°), and fields with a large object that 
covered a substantial portion of the field. This left us with 
1494 fields. Because of repeat visits by the HST separated 
by more than 6 months and visits to fields that are sep- 
arated by only of order tens of arcseconds, some of these 
fields overlap. So as not to double count fields, we dis- 
carded one of any pair of fields that lie closer than 1' to 
each other (the STIS Field of view is 50" x 50"). When 
considering which field to discard, we kept the field with 
the greater total exposure time. 

Due to several reasons including thermal effects, the 
HST undergoes jitter, or shaking (Nclan & Makidon 
2002). This jitter tends to blur images and introduces 
changes into measured object shapes. The median root 
mean square (RMS) jitter in one direction is 3-4 mas, but 
the RMS jitter value can be greater than 100 mas. We 
thus discard fields in which any of the input exposures 
had an RMS jitter greater than 10 mas (0.2 STIS pixels) 
in either of two orthogonal directions, or for which no jit- 
ter information is available. After these cuts, we were left 
with 1292 fields. 

3. PROCEDURE 
3.1. Shape Measurement 
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We follow the method described in RRGI to measure 
the object shapes and sizes and to correct for systematic 
effects. We measure the second Iki and fourth. Ikimn 

order 

moments for each object. These moments are defined as a 
sum over pixels p given by, for example, 
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Y, p x k xil(x)w(x) 

E„/(xMx) 



(2) 



where /(x) is the intensity in a pixel with position x = 
(xi,x 2 ) — (x,y) with respect to the centroid of the ob- 
ject and w(x) is a Gaussian weighting function. Similar 
equations hold for the remaining moments. The standard 
deviation of the weight function is chosen to be (in pixels) 



a w = max 




(3) 



where the area A is determined from the SExtractor 
measured semi-major(a) and semi-minor(6) axes as A = 
7r ( £ ^) 2 - The minimum weight function size of 3 pixels 
(0.15 arcseconds) is the optimal weight function we found 
for point sources (stellar images). This is similar to the 
minimum weight function size of 2 pixels (0.2 arcseconds) 
found in RRGI for WFPC2. 

These measured moments are used to compute the two 
component ellipticity of each galaxy given by 
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The first component of this ellipticity (ei) indicates elon- 
gation along the x (positive ei) and the y (negative e\) 
axes, e-i indicates elongation along axes 45° and —45° 
from the x axis. We also define the RMS size d of the 
object as 

d 2 = l(hi+l22)- (5) 

There are several systematic effects which affect the 
measured shapes of galaxies. The PSF of the telescope 
introduces smearing that can be deconvolved into an 
isotropic and an anisotropic component. The optics intro- 
duce an instrumental shear into the galaxies. Our method 
to correct for these effects is presented in RRGI. Here, we 
describe the particular steps taken for the STIS data set. 

3.2. PSF Model 

Ground-based surveys typically use stars in the survey 
fields to measure and correct for the PSF. Due to the 
small field of view of HST, space-based surveys do not 
have enough stars per field to do such a correction and 
have relied on separate observations of high signal-to-noise 
(S/N) stellar fields to measure the PSF (see e.g. RRGI; 
RRGIII; Hammerle et al 2002). However, in examining 
the STIS PSF, we found that the PSF shape depends on 
the S/N. As shown in Figure 3, the average ellipticity (ei) 
of stars in the survey fields changes from negative (elon- 
gation along the y axis) to positive (elongation along the 
x axis) as S/N increases. This result is consistent with a 
previously known charge transfer efficiency (CTE) effect 
in STIS (Kimble, Goudfrooij, & Gilliland 2000). Objects 
with few counts (due to short exposure times or low S/N) 



bleed in the y direction causing an elongation of the object 
in that direction. This CTE effect has been shown to be 
worse in the (y < 512)half (hereafter the "bottom half") 
of the STIS CCD (Woodgate 2002). Figure 3 confirms 
that the elongation in the y direction is indeed more pro- 
nounced in the bottom half of the STIS CCD, providing 
further evidence that the CTE is indeed the source of the 
S/N dependence of the PSF. 
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Fig. 3. — Stellar ellipticity ei as a function of signal to noise 
(S/N) for stars in galaxy-dominated fields. The triangles are binned 
averages in the top half of the chip, squares are binned averages in 
the bottom half of the chip. The error bars represent 1 <r errors in 
the mean. The dashed lines represent linear fits for 20 < S/N < 100 
and quadratic fits for S/N < 20. These fits are used to derive 
corrections for the second order moments ij.; based on S/N and y 
position. 

The CTE effect described above complicates the mod- 
elling of the PSF, and thus the correction of galaxy shapes 
for PSF smearing, for two reasons. The first is that we can- 
not simply measure high S/N stellar PSFs and use these 
to correct low S/N galaxies. The second is that we cannot 
use separate stellar fields to measure the PSF. Stars in stel- 
lar dominated fields (such as images of globular clusters) 
do not suffer the same CTE effects as stars in sparsely 
populated fields, even at low S/N. This is because the 
higher background levels in stellar fields fill in the charge 
traps causing the CTE effect. Thus, low S/N stars in 
high-background fields have a PSF similar to the PSF of 
high S/N stars regardless of background. Since the bulk of 
our galaxies are low S/N galaxies in low-background fields, 
that is the regime in which we must model the PSF. 

For S/N> 100 we see no evidence for an e\ dependence 
on either S/N or y position. Therefore, we consider the 
high S/N PSF to be the "true" STIS PSF and model the 
S/N and y dependence of the PSF as perturbations about 
this true PSF. We select 534 high (> 100) S/N stars from 
our survey and fit the 8 stellar moments from these stars 
to 5th order polynomials in x and y position as described 
in RRGI. Since we are attempting only a first order correc- 
tion for the effects of S/N and y position, we do not correct 
the fourth order moments for these effects. The low S/N 
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moment at any position (x,y) is given by the high S/N 
moment at that position, as calculated from the 5th order 
polynomial fits, minus some correction factor. The cor- 
rection factor is the difference between the average value 
of the moment at high S /N and the average value of the 
moment at low S/N at the y position of the star on the 
chip. Thus, 



I l k T(x,y) = I h k l 9h (x lV ) - - I l k T(y, S/N)), (6) 

for each diagonal component k = 1,2. 
We then write the moments as 

In - (l + ei )rf 2 , I 2 2 = {l-e 1 )d 2 (7) 

where the size d is defined as in equation (5). This size 
is slightly dependent on S/N, but has no dependence on y 
position. The size d 2 as a function of S/N is well fit by a 
linear equation given in Table 1. We fit e\ to the S/N in 
two different bins of y position, y < 512 and y > 512. We 
use linear fits for the range 20 < S/N < 100 and quadratic 
fits for the range 8 < S/N < 20. Table 1 gives the best fit 
functions, which are plotted on Figure 3. To calculate the 
average value of e\ at a particular S/N and y, we perform 
a linear interpolation or extrapolation between the two y 
fits at that S/N value assuming that the y < 512 fit is valid 
for y — 256 and the y > 512 fit is valid for y = 768. We 
find that there is no S/N or y dependence for e2, which is 
consistent with zero. Thus, no 1 12 correction is necessary. 
These fits allow us to calculate the corrections to the high 
S/N stellar moments using Equation (6). 

The average stellar ellipticity we measure at high S/N 
(S/N> 100), 3-5%, is slightly higher than the 1-2% mea- 
sured by Hammerle et al (2002). We compare our mea- 
surement with theirs by obtaining one of the high S/N 
stellar fields (o48b41010_3_ass) in their publicly available 
catalog and performing our analysis on that image. This 
image is drizzled to increase resolution and thus has a lin- 
ear pixel size of 0.025" (half of the intrinsic STIS pixel 
size) . If we analyze this image using the method outlined 
above with a Gaussian weight function width a = 6 pixels 
(twice what we use in our non-drizzled images), we find 
an average stellar ellipticity e\ — 4 — 5%. This is stable 
for larger weight function widths, but if we use a smaller 
weight function width (a — 3 drizzled pixels) we measure 
an average e\ — 1.5%. Thus, we are able to recover the 
value measured by Hammerle et al (2002) by using a suffi- 
ciently compact weight function. The weight function we 
use to measure stars represents a trade-off between the 
reduction of noise in moment measurements and the mea- 
surement of the full PSF anisotropy. As in RRGI, stellar 
moments are relatively stable for small perturbations in 
this weight function width. 

3.3. Correction for the PSF 

Using the stellar moments-corrected for weighting and 
detector shear as described in RRGI and fit to S/N and y 
position as described above-we can correct the galaxy mo- 
ments for the effects of PSF. The first step is a correction 
for the anisotropic smearing of the PSF and the second 
step is a correction for the isotropic portion of the PSF 
which increases the galaxy size slightly. This correction 
method was specifically tuned to HST images with small 



PSF and has produced excellent results in the correction 
of WFPC2 images for PSF effects (RRGII; RRGIII). The 
corrected galaxy moments are used to calculate the ellip- 
ticity of each galaxy, which allows the calculation of the 
overall shear in each field as described below. 

3.4. Shear Measurement 

We detected 63,895 objects in the 1292 fields we used. 
We discard stars and objects for which no jitter informa- 
tion is available or for which the jitter is very high. We fur- 
ther discard objects which have spurious ellipticities (non- 
physical post-correction ellipticities usually caused by poor 
background subtraction), leaving us with 49692 objects. 
We then discard objects that are too small to allow a shape 
measurement, i.e. with d < 1.7 pixels (after corrections), 
leaving 35033 objects. In RRGII and RRGIII, a magni- 
tude cut was made to discard galaxies too faint to allow 
for shape measurement. Here, we opt to use a S/N cut 
based on SExtractor values of flux and flux error. A cutoff 
value of S/N = 10 left us with 26036 galaxies. 

We thus obtain the corrected ellipticity of each se- 
lected galaxy that passed all the cuts mentioned above. 
We then calculate the shear estimator 73 for each galaxy 
using the equation 



where the shear susceptibility factor (G) ~ 1.4 is calcu- 
lated according to Equation 28 of RRGI and averaged over 
all galaxies in our final sample. 

4. RESULTS 

4.1. Shear Measurement and Cosmological Parameter 
Estimation 

To quantify the cosmic shear statistics, we follow the 
procedure of RRG III. The cosmic shear variance cf cns is 
obtained from a weighted sum of the estimators <rf ens j for 
each field / as 

ofens = ^2w f af ensJ , (9) 
/ 

where the weights Wf oc cr~ ^ se ^ are normalized so that 
J2f w f = 1- We define the noise variance t7„ oise j = 
(|7j olse | 2 ), which is measured from the data by computing 
the error in the mean 7^ from the distribution of galaxy 
shears in each field (RRGIII) . This provides a nearly opti- 
mal weighting scheme which allows us to use all fields that 
are not contaminated by stars or bright objects. Fields 
with few objects that pass our selection criteria are appro- 
priately down-weighted. 

With our STIS data set, we find 

of ens = (5.43 ± 1.06 ± 1.74) x 10~ 4 , (10) 

where the first error corresponds to noise only, while the 
second includes noise, cosmic variance and systematics. 
Our detection of cosmic shear (first error) is thus signif- 
icant at the 5.1cr level. The systematics error is domi- 
nated by the uncertainty in our PSF correction (see §3). 
To estimate it, we first noted that the uncertainty in 
our measurement of the PSF ellipticity in figure 3 is 
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Ae* ~ 0.01. The resulting uncertainty in the galaxy shears 
is A 7 ~ G- l ((d*)/(d)) 2 Ae* ~ 1.0 x 10" 3 (see RRG I, 
equation 58), where the average rms size (see Eq. [5]) of 
the PSF (with S/N=10, see table 1) and of galaxies (af- 
ter our cuts) is (d*) ~ 1.2 pixels and (d) ~ 3.2 pixels, 
respectively. 

Our measurement of the shear variance is shown in 
figure 4. Recent measurements from other surveys with 
similar galaxy redshifts (z m ~ 0.8-1) are also shown for 
comparison (RRGIII; van Waerbeke et al. 2002; Brown et 
al. 2003; Bacon et al. 2002) are also shown for compar- 
ison. The solid lines show the predictions for a ACDM 
model with Q m = 0.3, erg = 1, and T = 0.21 and with 
the above range of galaxy median redshift z m . Our vari- 
ance measurement is consistent with these other measure- 
ments and with the ACDM prediction. It also agrees with 
the measurement of Hammerle et al. (2002) who found 
of ens ~ (15 ± 12) x 10~ 4 (from their figure 15). Note that 
their error bar is considerably larger than ours, because 
their measurement is based on about 1/10 of the fields we 
have used. 
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Fig. 4. — Cosmic shear variance a? as a function of the radius 9 

lens 

of a circular aperture. Our observed value (HST-STIS) is shown as 
well as that from other recent measurements: van Waerbeke et al. 
(2001), Bacon et al. (2002, WHT and Keck), Brown et al. (2003), 
Refregier et al. (2002, RRGIII, HST-WFPC2). Also displayed are 
the predictions for a ACDM model with Q m = 0.3, T = 0.21, <rg = 
0.7 and 1.0, and a range of galaxy median redshift z m = 0.8, 0.9 and 
1.0, using the Smith et al. (2003) nonlinear fitting function. When 
relevant, the inner error bars include statistical errors only, while 
the outer error bars include statistical errors and cosmic variance. 



Our measurement allows us to set constraints on as, the 
amplitude of matter fluctuations on 8/i -1 Mpc scales. In 
a ACDM model with n m = 0.3 and T = 0.21 the shear 
variance is given by 



lens 



5.17xl0" 4 -4 



^8 
1.0 



£7 m 
03 



1.10 



0.21 



0.44 



\i.oJ 

where z m is the median redshift of the galaxies, whose 



with zq — z. m /1.09. The scale 0.51' is the effective radius 
of a circular cell corresponding to the chip size (50" on a 
side). The details of this calculation can be found in Bacon 
et al. (2000). Note however, that unlike these authors, we 
have used the more accurate fitting function of Smith et 
al. (2003) to compute the non-linear corrections to the 
power spectrum. 

For our sample and weighting scheme, the effective me- 
dian magnitude of the galaxies is M m ~ 24.8. According 
to Equation 1, this corresponds to a median redshift of 
z m = 1.0 ± 0.1, where the la error arises from the un- 
certainty in this equation. Our measurement of of ons thus 
yields 



cr 8 = (1.02 ±0.16) 



0.3 



0.46 



0.21 



(12) 



redshift distribution was assumed to be p(z) 



~(z/ z o) 



where the error includes noise, cosmic variance, system- 
atics, and the redshift uncertainty. This result is consis- 
tent at the lcr level with some earlier cosmic shear sur- 
veys (Bacon et al. 2002; RRGIII; van Waerbeke et al. 
2002) which derived values of er 8 in the range 0.9-1.0 for 
fi m = 0.3 (see figure 9 and table 1 in Refregier 2003 for 
a summary of the latest cosmic shear results). It is also 
consistent, but only at the < 2a level, with other cos- 
mic shear surveys which found erg in the range 0.7-0.8 
for the same value of f2 m (Brown et al. 2003; Hamana 
et al. 2003; Jarvis et al. 2002; Hoekstra et al. 2002). 
The recent WMAP results (Spergel et al. 2003) yield 
a 8 = 0.91 ± 0.21(0.3/^ m ) 6 when taken alone (for their 
single power law model), and a s = 0.78±o °o(0.3/rj m )°- 6 
when combined with other CMB and large-scale structure 
data (for their running power law index model). Our re- 
sults are thus also consistent with these WMAP values at 
the < 1.5ct level. 

4.2. Temporal Stability of the Shear Measurement 

The CTE effect described in section §3.2 is growing 
worse as the STIS CCD ages (Goodfrooij & Kimble 2002). 
In order to test whether this degradation affects our re- 
sults, we perform two tests. The first test involves stars 
and the second test involves examining our final result for 
temporal stability. 

There are sufficient stars in the survey to perform the 
analysis described in §3.2 on two subsections of the data. 
We divide the data roughly in half by time, classifying data 
taken before mid 1999 as "early" and data taken after that 
time as "late." We find that for S/N> 20 there is little 
change to the fit to e\ shown in Figure 3. However, for 
S/N< 20 we find that the fits shown in Figure 3 do change 
slightly. The early data are fit by curves with e\ slightly 
higher than those shown in Figure 3 and the late data show 
a slightly lower e±. This is consistent with a CTE that is 
growing worse over time. In both cases, the departure from 
the curve shown in Figure 3 is about 1%. This average 1% 
systematic error in the correction of stellar ellipticities at 
low S/N is included in the error estimates for our values 
of a^ and erg (see §4.1). 

To perform the other test, we divide our galaxy data 
into nine equally spaced date bins. We then measure the 
RMS shear a^ in each bin. We show these measurements 
in Figure 5. This plot shows no trend over time as might 
be expected if the CTE effect was not being sufficiently 
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corrected for. The shear signal is consistent throughout 
the time in which the data we use in this paper was taken. 
The only marginal outlier is the very first bin, which is 1.7cr 
from the value for the entire survey when statistical errors 
and cosmic variance are included. We have run several 
tests to determine why data taken during this period gives 
a higher shear signal than data taken during the rest of the 
survey but find nothing special about the data in this bin. 
Our results in this date range are consistent with those 
of Hammcrle et al. (2002) who found a high shear signal 
using STIS data from roughly the same period. Excluding 
the first bin (by only analyzing data taken after December 
1997) only changes the measured a^ from 5.43±1.74x 10 -4 
to 4.76 ± 1.68 x 10" 4 . The effect on £78 is less pronounced; 
the removal of the first bin decreases as from 1.02 ± 0.16 
to 0.97 ±0.16. 
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Fig. 5. — Cosmic shear variance o"j 2 ens as a function of the date. 
The inner error bars include statistical errors only, while the outer 
error bars include statistical errors and cosmic variance. The dotted 
horizontal lines represent the value of Cj 2 ens and the associated one 
sigma statistical error for the entire survey. The left-most point is 
consistent with what Hammcrle et al. (2002) found using roughly 
the same data. 



5. CONCLUSIONS 

We have measured the shapes of 26036 galaxies in 1292 
STIS fields, corresponding to about 0.25 square degrees. 
We corrected for systematic effects as outlined in RRGI 
with the added step of correcting galaxy moments for a 
S/N dependent charge transfer efficiency. We have de- 
tected a cosmic shear signal using the STIS parallel archive 
at the 5.1(7 level. After correcting the galaxy shapes for 
PSF distortions, detector shear and CTE effects, we mea- 
sure a shear variance of <r\ cns = (5.43 ± 1.74) x 10~ 4 on 
the STIS chip scale (50" on a side), where the la er- 
ror includes noise, cosmic variance and systematics. This 
is consistent with the earlier and much noisier measure- 
ment of Hammcrle et al. (2002). For a ACDM model 
with fl m = 0.3 and T = 0.21, this sets a constraint on 
the amplitude of the matter power spectrum given by 
a s = 1.02 ±0.16, where the la error includes noise, cosmic 



variance, systematics and redshift error. This is consis- 
tent with earlier measurements of as from cosmic shear, 
but tends to favor those with higher values. It is also con- 
sistent with the recent determination of as (Spergel et al. 
2003) from CMB anisotropies with the WMAP mission. 

The results presented here represent a contribution to 
the first generation of space-based weak lensing results pre- 
sented in RRGII, RRGIII, & Hammerle et al. (2002). The 
next generation of space-based parallel observations opti- 
mized for weak lensing will be made with the Advanced 
Camera for Surveys (ACS) on the HST. A targeted ACS 
project, the COSMOS 2 square degree field (GO-9822) will 
provide strong constraints on cosmological parameters and 
allow for high resolution dark matter map. ACS surveys 
will benefit from the high resolution of ACS (0.05 arcsec- 
ond pixels like STIS) as well as its enlarged area (approx- 
imately 10 square arcminutes) and improved sensitivity 
(relative to WFPC2). New methods for the measurement 
of object shapes and the correction of systematic effects 
are being developed that will capitalize on the excellent 
resolution of this survey and future space-based surveys 
(Refregier & Bacon 2003; Bernstein & Jarvis 2002). Fu- 
ture generations of weak lensing surveys both from the 
ground (e.g., the Canada- France-Hawaii Legacy survey, 
Mellier et al. 2000; the Large-aperture Synoptic Survey 
Telescope; Tyson et al. 2003) and from space (e.g. the 
Supernova Acceleration Probe; Alcock et al. 2003; Rhodes 
et al. 2003, Massey et al. 2003; Refregier et al. 2003) will 
continue to utilize this unique method to measure cos- 
mological parameters to unprecedented accuracy. These 
projects will survey large areas while obtaining images in 
multiple filters allowing for accurate photometric redshifts 
of the survey galaxies. SNAP, with its near-infrared capa- 
bility, high resolution, and wide field of view, will be able 
to probe mass concentrations beyond z = 1 , thus allowing 
for redshift tomography and the study of the growth of 
structure. This will enable weak lensing to make a contri- 
bution to the study of dark matter and dark energy, and 
to set tight constraints on cosmological parameters. 
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Table 1 

PSF Model Parameters for Equations 4 and 5 



Parameter 


S/N range 


y range 


Fit 


d 1 


10-100 


all 


1.47 + 0.00072 x S/N 


ei 


8-20 


y < 512 


0.124 - 0.023 x S/N + 0.00089 x (S/N) 2 


ei 


8-20 


y > 512 


0.010 - 0.0025 x S/N + 0.00021 x (S/N) 2 


ei 


20-100 


2/ < 512 


0.016 + 0.0022 x S/N 


ei 


20-100 


y > 512 


0.031 + 0.00015 x S/N 


£2 


all 


all 






